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The Schrodinger equation for an interacting spinless electron gas in a nonuniform magnetic field 
admits an exact solution in Jastrow product form when the fluctuations in the magnetic field track 
the fluctuations in the scalar potential. For tracking realizations in a two-dimensional electron gas, 
the degeneracy of the lowest Landau level persists, and the "tracking" solutions span the ground 
state subspace. In the context of the fractional quantum Hall problem, the Laughlin wave function 
is shown to be a tracking solution. Tracking solutions for screened Coulomb interactions are also 
constructed. The resulting wavefunctions are proposed as variational wave functions with potentially 
lower energy in the case of non-negligible Landau level mixing than the Laughlin function. 



I. INTRODUCTION 

Modern semiconductor technology has produced ex- 
ceedingly pure materials, and nanostructures that have 
a number of interesting properties at low temperatures. 
A comprehensive review is given by Beenakker and van 
Houten.El Modulation-doped GaAs-Al^Gai-^As hetero- 
junctions with precisely defined narrow potential wells at 
the interface strongly quantize electronic motion pepen- 
dicular to the interface, creating a two dimensional elec- 
tron gas (2-DEG). If a magnetic field B a is applied per- 
pendicular to the interface, electronic motion is further 
quantized in the interface plane, giving rise to a series 
of condensed Landau energy levels. At high magnetic 
field values the Hall conductance is accurately and ro- 
bustly quantized in integer multiples of e 2 /h. The re- 
sulting "integer quantum Hall effect" can be understood 
to a large extent within a one-electron picture. On the 
other hand, quantization of the Hall conductance at frac- 
tional filling factors, the "fractional quantum Hall effect" 
(FQHE), is intrinsically a many-body effect in which 
electron-electron interactions play the principal role. 

Near Landau filling factors v — n/(2mn + 1), with 
n, m integer, the Hamiltonian of the strongly interact- 
ing 2-DEG can be transformed to a weakly interact- 
ing composite-fermion (CF) model in which each elec- 
tron is coupled to an even number of magnetic flus 
quanta </> = 2m. Since its introduction, the CF.-mpdelEl 
has attracted considerable theoretical attention.ElLl The 
essence of the electron-electron correlations is captured 
by the single-particle mean field approximation for the 
CF's, which are subject to a reduced magnetic flux den- 
sity B c s = B a (l — <j)v). Since the CFs observe a re- 
duced mean magnetic field, they fill up a modified Lan- 
dau level structure, in which the degeneracy of each level 
is reduced appropriately. As a result, the weakly inter- 
acting CF gas is described by an effective filling factor 
z'cff — vB a /B D f[ which becomes integer z/ c ff = n at the 
above special fractional v values. The fractional quantum 



Hall effect follows by assuming that the CF's manifest the 
integer quantum Hall effect when i/ g is integer. 

The composite fermion (and similar composite bo- 
son) theories of the fractional quantum Hall effect were 
largely motivated .by the success of Laughlin's variational 
wavefunction.EJO Laughlin's wavefunction attaches col- 
lective vorticity to the electron gas which efficiently re- 
duces the probability of two electrons approaching each 
other, and the Coulomb energy is reduced when v is a 
simple odd fraction v = l/(2m + 1). It was soon realised 
that such vorticity can be generated by a singular gauge 
transformation which attaches infinitesimal flux-tubes to 
the particles, i.e., transforms electrons into composite ob- 
jects with transmuted statistics.u However, such singular 
gauge transformations by themselves are not capable of 
reproducing the full Laughlin wavefunction, since they do 
not give its correct radial dependence. The full Laughlin 
wavefunction is obtained by taking into account fluctua- 
tion effects around the mean field solution.tj 

The Laughlin wavefunction is a trial state completely 
formed within the subspace of the lowest Landau level. 
It is believed to reflect the most important aspects of the 
FQHE ground state at odd fractions, v — l/(2m+ 1), in 
the high magnetic field limit, where Landau level mixing 
can be neglected. Similarly, suggested generalizations to 
other fractional filling factors such as v = n / (2mn+l) are 
appropriately projected onto the lowest Landau level.El 
In the ground state it is a priori not clear how to take 
a strong electron-electron interaction into consideration, 
or alternatively a weaker magnetic field, when the mixing 
of the Landau levels is not negligible. In such cases the 
true many-body ground state also has components from 
higher Landau levels. 

In this paper we develop exact eigenfunctions \& of the 
Schrodinger equation for an interacting electron gas in 
non-uniform magnetic fields. The initial insight and mo- 
tivation is taken fcpm the exactly s&Lvable ID construc- 
tions of CaIegerot3 and Sutherlandca, PolychronakoalJ, 
and Halcytj in the absence of magnetic fields, and from 
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the 2D constructions for the Eauli Hamiltonian intro- 
duced by Aharonov and GaaheiO, and further developed 
by Dubrovin and NovikovO for an electron in a periodic 
magnetic field. Our analysis is based on an ansatz that 
linearizes the vector Riccati equation in S = — Vln^, 
which is equivalent to the Schrodinger equation. Our 
linearization is valid in any number of dimensions, but it 
is most easily utilized in two dimensions. We refer to the 
ansatz based solutions as "tracking" solutions, since they 
satisfy certain linear first order differential equations con- 
taining the vector potential, subject to a constraint re- 
lating the vector and scalar potentials. We demonstrate 
that the Laughlin wave function and any other wave- 
function completely in the subspace of the lowest Landau 
level, is a tracking solution. This tracking solution uti- 
lizes the Chern-Simons vector potential a cs , but in a dif- 
ferent way than used in the composite fermion theories: 
In our theory, a cs does not appear in the Hamiltonian. 

We also derive a new trial wave function that is 
the exact tracking solution for the screened Coulomb po- 
tential. A priori, is just as close to the real physical 
situation as ^l- While misses the correct treatment 
of the Coulomb interaction, ^fc requires the presence of a 
two-body magnetic field that tracks the scalar potential. 
The advantage of $c, however, is that when it is used 
as a trial wavefunction, it depends on free variational pa- 
rameters and thus potentially gives a lower energy than 
^ l, which has no adjustable parameters. The wavefunc- 
tion ^>c is proposed to improve the FQHE ground state, 
especially in the case when Landau level mixing is im- 
portant. 



II. GENERAL FORMULATION AND TRACKING 

Consider the time-dependent Schrodinger equation for 
N electrons, each with charge — e < and effective mass 
m*. An electron with coordinate is subject to a scalar 
potential V(r.i,t) and a vector potential A(rj,t), which 
in general include contributions from both external fields 
and internal interactions depending on the coordinates of 
all the electrons. In this context, the Schrodinger equa- 
tion has the form 
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(2.1) 



where {r} = {ri, ra, ...rj, ...rjv}, and Vj is the gradient 
with respect to the electron coordinate iv The constant 
7 = 2m* /H 2 , and 0o = h/eis the flux quantum. Particles 
are labeled by Latin subscripts, but the constant i = 
v 7 — 1. To assure the hermiticity of the Hamiltonian, the 
vector potential A(r i; t) is assumed to be real. The scalar 
potential has the form 



V(r«, t) = V a (Ti,t) + V c (ri - 



(2.2) 



where V a [vi , t) is a time-dependent applied scalar poten- 
tial, and V c is an internal two-body electron-electron in- 
teraction potential, such as the Coulomb interaction. 



A. Riccati Equation 

The starting point of the ana lysis is a reformulation 
of the Schrodinger equation (2.1) in terms of a function 
Q({r},t) defined by 



*({r} J *)=exp[-Q({r},t)]. 



(2.3) 



In general, Q = — In $ is a complex multivalued function. 
Substituting W from Eq. (2.3) into Eq. (2.1) leads to the 



nonlinear vector Riccati equation 
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where 



S 4 = S 4 ({r},i)=ViQ({r},t), 



(2.5) 



Aj = A(r.;,i), and Vi — V(ri,t). The above form of 
the Riccati equation contains $ as a multiplicative fac- 
tor on both sides. In the following this factor will be 
divided out, even though such division is not legitimate 
at the zeros of VP. This subtlety results in the appear- 
ance of spurious delta-function potentials at the position 
of the zeros. However, as shown in section II, the correct 
physics is preserved. 

Although the original Schrodinger equation is a linear, 
second order equation for \E', Eq. (^J) is a first order, 
nonlinear equation for the derivati ves of Q. Given real 
vectors Aj and real scalars Vi, Eq. (|2.4| ) can, in principle, 
be solved for Q({r},t). For time-independent potent ials , 
the function Q({r}, t) = Q({r}) + i(E/h)t, and Eq. ( jgj ) 
becomes the time-independent energy eigenvalue equa- 
tion 
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In the following we restrict our attention to time- 
independent problems, but it will be noted that the track- 
ing constru ction introduced in subsection II. B also lin- 
earizes Eq. (2.4). 

In one d ime nsion the transformation defined by Eq. 
J3h and (2.1) leads to a scalar Riccati equation that 
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is quite useful for obtaining the stationary, one-particle 
states for many scalar potentials in the absence of mag- 
netic fields.li^l In the ma ny-particle case the major diffi- 
culty in solving Eq. (2.6) is that the quadratic term in Si 
may contain interactions of order higher than those ap- 
pearing in the scalar potential Vi- To see this, consider 
wave functions in a Jastrow product form, 



*(w)=nn n ( r *- r i)-n^ 



(2.7) 



where Q(rj — rj) is a two-body function depending on the 
coordinate difference between the particles i and j, and 
Xi = x( r i) is a one-body function involving only the co- 
ordinates of particle i. Due to the two-body factors fi^, 
the vector Sj will be a sum over two-body functions, each 
of which depends on the coordinates of the ith electron 
and another one; thus the quadratic term in Si implicitly 
contains three body interactions. 

One may conclude that a wavefunction in the Jastrow 
product form of Eq. (ET7|) can only be an eigenstate if 
the model in question contains an unphysical three-body 
potential. Boson Hamiltonians of this kind have been in- 
vestigated in Ref. There it is argued that even though 
the three-body term is explicitly present in the Hamilto- 
nian, it becomes irrelevant in the renormalization group 
sense; thus the Jastrow product, although it is not an ex- 
act eigenstate, describes correctly the low-energy, long- 
wavelength physics of the remaining two-body Hamilto- 
nian. While three-body terms are required in the generic 
case, Jastrow products can be exact eigenstates of Hamil- 
tonians without three-body terms, provided that special 
conditions are fulfilled, as shown in section II. B. 

Well-known examples in—piic— dimension are the 
Calogero-Sutherland models .E3i1j'l3 They form a class of 
integrable many-body systems, which describe particles 
interacting through two-body potentials V c — r.j) in 
zero magnetic field A; = 0. Their g rou nd state can be 
written in the product form of Eq. fl2.7p. For such sys- 
tems the sum of the quadratic terms of the associated ID 
Riccati equation reduces to a sum of two body interac- 
tions 
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The function = F(n — Tj) is symmetric, and only de- 
pends on the coordinate difference between the ith and 
jth particles. This equation is a consequence of a more^- 
ementary functional equation (Sutherland's equationta) 
assuring that the scattering of three particles can be fac- 
torized into two -par ticle scatterings. The class of models 
satisfying Eq. (2.S) is however limited: Typically, the 
two-body potential V c is constrained to diminish over 
short distances as 1/r?-. 

The Riccati equation, Eq. (2.6), provides an alternative 
description of the quantum mechanical system. However, 
care should be taken since not every solution of the Ric- 
cati equation yields a valid, single valued quantum me- 
chanical wavefunction (we do not consider anyons here). 



Assume that a solution Si of Eq. ( |2.6D is found. The 
function Q can then be obtained by inverting Eq. ([2.5|), 
which is tantamount to the line integral 

,{*} 

Q({r}) = Q({r })+ {rfr'} • {S({r'})}, (2.9) 

J {ro} 

where {S} = {Si, S2, • • • , Sjv}, and {dr} = 
{dri, dx2, . . . , drjv} are DiV-component supervectors (D 
is the dimension of the space), and Q({i"o}) is a constant. 
In the many-particle case the line integral is path depen- 
dent whenever the DN dimensional vector field {S} is 
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non-conservative, i.e. 

i, j = 1, . . . , N; a, (3 — 1, . . . , D. Since the wavefunction 
must be single valued, the function Q is determined 
only up to an additive term z27tto, with m integer. Hence, 
we allow a set of well-defined, isolated, singularities of 5- 
function type in any two-dimensional slice (ia,j(3) of the 
DN dimensional space. Restricting our attention only 
to solutions in the Jastrow-product form, it is enough to 
require a constraint on V.; x Sj, i.e., for slices in the same 
particle subspace with i = j, since then the special form 



in Eq. (2.7) assures that, at most, similar <5-function type 
singularities appear in other slices with i =/= j. 

In 3D we allow discrete loops described by vectors 
Cfc(s), parametrized by the length s, on which x Si is 
proportional to a (^-function. We write 

V,xS i = -i2tt ^m fe % - c fc (s)) d s c k {s), (2.10) 

k 

where is an integer, and d s Ck(s) is a unit vector tan- 
gential to loop k at the point s. Although the singularity 
loops are independent of r^, they may depend on th e co- 
ordinates of the other particles. By virtue of (2.1C), the 
value of the line integral in Eq. (2.9) is determined only 



up to the term «27r^ fe rrikrik, where the integer rik is the 
Gauss linking number (winding number in 2D) of the 
path linked by the fc th singularity loop. Since J^ fc iriknk 
is int ege r, Eq. ( 2.10 ) assures that the wavefunction ^ in 
Eq. (2.3) is single valued. In a 2D plane defined by the 
perpendicular unit vector z, the vectors Cfe(s) — > c/. are 
discrete points, and d s Ck(s) — > z. 



B. Tracking Constructions 

In this subsection we construct a class of exactly 
solvable many-body Hamiltonians in higher dimensions, 
which contain both scalar and vector potentials. The 
many-bod y ei genfunctions can be written in the prod- 
uct form (2.7), with each term related to terms of an 
appropriate decomposition of the vector potential. Our 
construction is motivated by the need to linearize the 
Schrodinger equation to eliminate the three-body inter- 
action terms generated by the quadratic terms in and 
A,; in (2.6). Linearization may be achieved by requiring 
that 
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S(S.-^A,) 2 = 0, 



(2.11) 



which is a sufficient condition for the factorizability of 
the wavefunction. Eq uatio n (2.11) is a stronger restric- 





tion than that in Eq. (2.8), since it requires Fij = 0. In 
ID with Aj = 0, realizations of Eq. (2.11) are extremely 
limited since this equation implies that the wavefunction 
is a coordinate independent constant. It is the presence 
of the vector potential, and especiall y the higher dimen- 
sionality which make the const raint ( 2.11 ) plausible. We 
show in the following that Eq. ( |2.11 ) provides nontrivial, 
useful solutions to higher dimensional models containing 
vector potentials. 



To manifest the implications of Eq. (2.11) we decom- 
pose the complex vector Sj as Sj = 3?Sj + i ^Sj a nd write 
out the real and imaginary parts of Eq. ( 2.11 ), respec- 
tively as 



V (5RS,) 2 - (3S, - ^A 
J^SiSi-^Si-^Ai) 



= 0, 



(2.12) 



For a system of identical particles, each term in Eq. ( 2.12 ) 
is zero, and consequently the two real vectors SftSj and 
3Sj — 27r/0oAi are orthogonal to each other and have 
the same length at any point of the configuration space 
{r}. Evidently, if only Aj is given, this restriction still 
allows an infinite number of possible solutions for the 
complex vector Si. To parametrize these solutions we 
formally decompose Aj into the sum of two components 
Aj = A] + A] 1 , with each ter m real . Equation (2.11), or 
equivalently the equations in (2.12), are satisfied by the 
ansatz 
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A(n) = \Ax A]\, 



(2.13) 



where n is a unit vector with arbitrar y dir ection, and 
the unit vector A] = A\/A]. Equation ( 2.13 ) completely 
specifies the vector Si in ter ms of n and the decomposi- 
tion of A. Using the ansatz ( [2.13 ), the Riccati Eq. (157" 
reduces to 
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where fJ-B = eh/2m* is the Bohr magneton. In general, 
the first term in Eq. (2.14) is a complicated nonlinear 
function of the angle between n and A 1 . Considerable 
simplification is achieved when n and the coordinate sys- 
tem can be chosen such that Z?j(n) and n are constants. 
In 3D these conditions are restrictive, but in 2D they 
can always be satisfied (See Section II). Assuming that 



Di{p) and n are constant, Eq. (2.14) reduces to the linear 
equations 
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and 
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E, (2.15) 



(2.16) 



Since the rhs of Eq.( 2.15) is a coordinate independent 
constant E, Eq. (2.15) can only hold if the spatial vari- 
ations of the component of the magnetic field B : (rj) = 
Vi x A] in the directi on n t rack t hose o f the scalar poten- 
tial. Although Eqs. (2.15) and (2.16) give a restriction 
only on the sum of N terms, in case of identical particles 
they hold identically for each particle i, and each vector 
potential A\ is in the Coulomb gauge, Vj • A\ = 0. 



By virtue of the single valuedness condition Eq. ( 2.10 ), 
whose rhs is purely imaginary, it follows using Eq. (2.13) 
that B* 1 = Vj x A] 1 is necessarily restricted to the form 
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(2.17) 



with m,fc an integer. Thus the construction (2.13) ex- 
cludes regular magnetic fields in B* 1 , only allowing sin- 
gular flux tubes proportional to (^-functions. It is evident 
that if we already possess a solution of the Schrodinger 
equation, we can construct a new Hamiltonian, with the 
associated new eigenfunction, by adding a term 5 A] 1 to 
the vector potential and the associated term <5QS; to 
the logarithmic deriva tive o f the wavefunction, providing 
that they satisfy Eq. ( 2.13|). T he reader recognizes that 



the imaginary part of Eq. (|2 .13) is related to gauge trans 



formations generated by real functions A(r,), which mod- 
ify the vector potential as SAi = VAj and th e pha se of 
t he w avefunction as £(e/7i)Aj. Note that Eqs. ( 2.13 ) and 
( p. 17 ) also allow singular gauge transformations, gener- 
ated by multivalued functions Aj. Singular gauge trans- 
formations may lead to more substantial modifications in 
the model than regular transformations by adding singu- 
lar flux tubes to the system and thus eventually trans- 
muting particle statistics. 

The single valuedness condition Eq. ( 2.1C| ) also implies 
a rest riction on the type-I component A+. Using Eq. 
( |2.13 ) in Vj x SRSj = 0, with Dj(n) and n constant, and 
Vj • A] = 0, gives the constraint 



(n-Vj)A? = 0. 



(2.18) 



Thus A] must be constant in the n direction. However, 
since Vj x A] is not restricted, regular magnetic fields 
are allowed i n B|. 

Equation ( 2.1 5| ) is the main result of our paper. It 
completely determines the decomposition of the magnetic 
flux density, and if a physical problem, defined by giving 
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the functions Vj or B*, and n admits tracking solutions 
then there is only one tracking sol ution up to possible 
gauge transformations. Equ ati on ( 2.15) together with 
the constraints in Eqs. fl2.16p - ) further im ply that 
the iV-particle Schrodinger equation in Eq. (2.1) has an 
exact tracking solution with energy E and eigenfunction 



e -K<2({r}) e -i3Q({r})_ 



(2.19) 



For identical particles, the problem reduces to a one- 
particle problem in external fields and in fields created 
by other particles. The TV-particle wavefunction can be 
written as a product of one-particle wavefunctions (or in 
some cases as a linear superposition of these) , where the 
one-particle functions also contain the coordinates of the 
other particles as parameters. 

The tr ackin g wave function ^> may be obtained by solv- 
ing Eq. (2.13), which expresses the real and imaginary 
parts of Sj = ViQ = — Vjln 1 !' in terms of the A] and 
Ap components, respectively, of the vector potential de- 
composition. However, \^\ can be expressed directly as 
a function of the scalar potential Vj: Assuming that the 
particles are identical and noting that A is real, Eqs. 
(|2.6P and (2.11) yield Poisson's equation 



(2.20) 



with the en ergy eigenvalue E determined by the tracking 
constraint (2.15). Since the tracking solution is a Jastrow 
product, we find that In |^| = — SRQ is determined by the 
sum of quadratures over the scalar potential, 

N . 

§RQ({r})= 7 £ J d °< GdM) [E/N-Vft)], 

(2.21) 



where Gu(ri,r^) is the J-di mens ional Green's function. 
Although we do not use Eq. (2.21) to obtain the tracking 
solutions in Section III, it is applicable and it may be 
useful in other problems. 

It is of interest to calculate the current density cor- 
responding to the tracking solution. The electron cur- 
rent density associated with some solution of the 
Schrodinger equation is defined as 
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(2.22) 



where Eq. (2.3) was used to obtain the latter form, and 
Vj = Sj — i^Aj is a normalized complex velocity. Ap- 
plying Eq. fl2.13p , the velocity is Vj = SRSi - z(27r/0 o )Aj, 
yielding the current density 



M = 



(2.23) 



Since the current is gauge invariant, the singular type-II 
component A] 1 does not contribute. 



III. TRACKING IN TWO DIMENSIONS 

In this section we restrict the interacting electron gas 
to two dimensions, and develop several examples of track- 
ing solutions. In 2D, with both Si and Aj restricted to 
the x—y plane, we seek tracking solutions with n = z. 
Noting that -Dj(z) = 1, in 2D the vector Si is defined by 
the equation 



= ViQ- — [-z 



A] + iA?] 



The tracking constraint (2.15) becomes 
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J2[»BB I (r i )+V(r i )]=E, 



(3.1) 



(3.2) 



with B\=i - (Vj x A\). The constraints ( |2.16| ) - ( |2.1S| ) 
are Vj • A\ = 0, 

Bf - z • (V, x Af) = -0o rrikSiri - c k ), (3.3) 



and dA\/dz = 0, respectively. 

In the following we assume that the particles are spin- 
less fermions, and that the system is subject to an applied 
uniform magnetic field B a = iB a , with B a > 0. The 
corresponding vector potential in the symmetric gauge, 
written in polar coordinates r = (r, (p), is 



-B a rip. 



(3.4) 



We now develop a few interesting solutions of the tracking 
relation Eq. (3.2), relevant to the fractional quantum Hall 
effect, which satisfy the constraints following Eq. ([3.2]), 
and construct the associated eigenfunctions. 



A. Non-Interacting Electrons 

When A = A a is the only potential which appears 
in the Hamiltonian, i.e., the scalar potential V = 0, the 
one-electron states form highly degenerate Landau levels. 
Restricting our attention to filling factors v — N<f>o /B a < 
1, we consider exclusively the lowest Landau level (LLL), 
which has the one-particle energy E/N — [iBB a . The 
single particle (N — 1) Schrodinger equation with vector 
potential A = A a has eigenfunctions with anti-analytic 
prefactors of the form 
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4> m (r)=c m Z m exp[ —X ZZ , Z = x+iy 



(3.5) 



is 



where l~ 2 = eB a /h = (2n/(j) Q )B a , and c m cx ;~( m + 3 / 2 ) 
a normalization constant. These one-particle basis func- 
tions ip m with definite canonical angular momentum mh 



5 



(to a nonnegative integer) span the LLL. Substitution of 
if} m into Eq. (2T) yields the corresponding energy eigen- 
value E m = fiBB a . Using ij} m and A a in Eq. ( 2.22 ) gives 
the circulating current density 



liLi\r,„\- ( ^ - l a '~r . 



(3.6) 



Since c 2 n oc Ba 1+3 ^ 2 , the current j m — ► in the limit 
B a — > 0, as expected when there is no force on the elec- 
tron. 

We now show that any basis state tp m is a tracking 
solution corresponding to a S function potential V(r, to). 
Using the formalism of section II with N = 1 , we decom- 
pose the vector potential A a as A a = A 1 + A 11 with 



A 1 = A a - a(r, to), A 11 = a(r, to), 



(3.7) 



where the particle index has been omitted, and a is a 
singular vector potential defined by 



It is of interest to point out that another tracking so- 
lution, associated with the decomposition 

A 1 = A a — a(r, to), A 11 = — a(r,m), (3.11) 

yields the same LLL ener gy eigenvalue E m = fisBa and 
current density Eq. (3.6), but with basis function tp m 
whose prefactor is analytic. The form ip m , which is the 
one used by Laughlin,til corresponds to a vector poten- 
tial A = A a — 2a(r, to). Hence there is a (negative) net 
flux tube attached to the origin of the coordinate sys- 
tem, which the moving electron observes and picks up an 
Aharonov-Bohm phase. 

Now co nsid er the TV-electron problem. The con- 
struction (3.7) can be generalized, using the particle- 



dependent decomposition 

A] = A -a(ri,TOi), A] 1 = a(r i; mi), (3.12) 



with a defined as in Eq. (3.8). We find that Q = 
1/4Z~ 2 J^i ZiZi — J^i m i lnZj, which leads to the eigen- 
function 
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The auxiliary vector potential a(r, to) attaches a ficti- 
tious flux tube of density 



b(r, to) = z • V x a = m<j) 5(r), 



(3.9) 



to the origin of the coordinate system, but since the flux 
tubes are attached with opposite sign in A 1 and A 11 , they 
cancel out in the Hamiltonian and have no direct physical 
effect. 

To form the corresponding tracking eigenfunction we 
solve Eq. flOj), with A 1 and A 11 given by Eq. ( pT^ . The 
imaginary part gives = mip, while the real part yields 
SRQ = l/4 Z~ 2 r 2 — to In r. Adding these results and using 
Eq. ( 2.19| ) gives the wavefunction ip m of Eq. (3.5). The 
current de nsity of the tracking sol ution can be calculated 
using Eq. ( 2.25 ), which yields Eq. (3.6). In the analysis it 
is seen that the A 1 piece is responsible for the the radial 
part of the polynomial prefactor in ip m) while the A 11 
piece gives the angular part. 

In order to satisfy the tracking equation Eq. ( |3.2| ) with 
energy E = /iBB a , we must introduce the scalar potential 



V(r, to) = Msfr(r, rn) = m^ B cj) 8(v). 



(3.10) 



Since the wavefunction ip m has a zero at the origin for 
to > 0, such a singular potential is absolutely harm- 
less and it does not contribute to the physical properties 
of the system. The appear ance of the delta function is 
due to the transformation (2.3), which transforms alge- 
braic zeros of ^ into logarithmic sin gularities of Q, and 
the subsequent division of Eq. (2.4) by the wavefunc- 
tion to obtain the time independent Riccati equation Eq. 
(2.6). Since V 2 -0 m is regular, the tracking solutions found 
also satisfy the original Schrodinger equation without the 
above <5-function potential. 



In the case of a fully filled LLL, v — 1, for a cir- 
cular droplet with unit area, the particle number is 
N = B a /<f>o, and the one-electron states are filled with 
increasing to quantum numbers up to N — 1. The 
total canonical angular momentum M = Ylj—-\ m i = 
+ 1 + . . . + (N - 1) = N(N - l)/2 is a conserved quan- 
tity. The angular momentum, however, can be arbitrarily 
distributed among the particles; thus any linear superpo- 
sition 



*M = ^ C{m}^{m}i 
{rn} 



(3.14) 



subject to the constraint J2i m i — -W) is a valid solu- 
tion. The A-electron wavefunction must be antisymmet- 
ric with respect to the exchange of the particles. This 
requires that C{ m } be the totally antisymmetric tensor 
of N indices, i.e., we have to form a Slater determinant. 
A closed forpi is achieved by forming the Vandermonde 
determinantE3 



{m} 



'{rn} 



det 



Zi 

Z 2 



1 z 



N 



21 



7 2 
N 



N-l 



z. 



N-l 



yN-l 
J N 



N N 

=n n^--^)- (3-i5) 

i=l j>i 

Thus the TV-electron wavefunction for the fully filled LLL 
takes the form 

N N / \ 

* m = n n & - z ^ cx p —fi 2 e z ^ > ( 3 - 16 ) 

i=l j>i \ i / 

where 2M = N(N -1). 
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B. Interacting Particles - The Laughlin Wave 
Function 

The above discussion remains equally valid if the sys- 
tem is subject to internal helds created by the particles 
themselves. Our general technique to find solutions to 
the many body problem is to solve the one-body prob- 
lem, form product wave functions from these solutions, 
and finally antisymmetrize them for fermions (or sym- 
metrize for bosons). 

Consider the addition of a Chern-Simons vector po- 
tential, whose value at the position of the ith electron is 
determined by the positions of the other particles as 



N 



z x 



(3.17) 



with cf> = (f> cs /(cj)QN) the number of flux quanta per parti- 
cle arising from a cs . The Chern-Simons vector potential 
a cs attaches fictitious flux tubes of density 



b cs (r,v 



N 



Vj x a cs (ri) = z(j)(j)Q Siji - r, 



(3.18) 



to the particles. These composite particles remain 
fermions if <p is even, while their statistics transmute to 
bosons when <f> is odd. 

As above, we decompose the vector potential A as 



A] = A a - a(rj, m) - a cs (r* , . 
A" = a(rj , to) + a cs (r, , cj>) . 



(3.19) 



Again the flux tubes appear with opposite sign in A 1 
and A 11 , so they cancel out from the Hamiltonian. The 
associated eigen f unction is form ed by solving Eq. (3.1) 
using Eq. (3. IE). Equation ( |3.2[ ) is satisfied with energy 
E = N^bBo,, provided that we introduce the fictitious 
scalar potential 



V(r,-,m, 



N 



)]■ (3-20) 



Recalling the discussion after Eq. (3.10), this singular 
potential is spurious and has no real effect on the prob- 
lem. Again, the wavefunctions to be found below are 
valid solutions of the original Schrodinger equation with- 
out V{vi,m, (f>). 

The Chern- Simmons contribution gives 



tan 



i j=jLi 



-i ( Vi - Vj 



3h 



which, combined with the other contributions calculated 
in the free particle case, gives 



Follow ing the anti-symmetrization procedure in Eq. 
( |3.15 ), the wavefunction takes the form 



i j>i 



3>i 



= n n & ex p ( - ?- 2 e . 

(3.21) 



where ^ m is the free electron wave function of Eq. (3.16). 
This is the Laughlin wavefunction, used to describe the 
state of the electrons for fractional filling factors 1 / (0+1) 
(with <p even) . The distinction between type-I and type- 
II potentials is that the type-I piece generates the radial 
part of the Laughlin prefactor, while the type-II piece 
generates the angular part. It is the type-I vector po- 
tential in the above construction that keeps the particles 
further away from each other, and increases the effective 
filling factor to v c fi = 1. 

Our derivation of the Laughlin wavefunction deserves 
some interpretation. It is well-known that without the re- 
pulsive interaction of the particles any wavefunction with, 
an anti-analytic prefa ctor is an eigenstate in the LLL.E2I 
The state in Eq. ( 3.21 ) is only one of the many degenerate 
ground states. Switching on the weak Coulomb interac- 
tion slightly mixes up the different Landau levels. The 
Laughlin function ^ l ceases to be an exact ground state 
but in the high field limit it is believed to become a good 
variational ground state, better, in general, than other 
wavefunctions with anti-analytic prefactors. 

Usually, the Laughlin wavefunction is derived in a field 
theoretic framework beyondj-thc mean ficM level. Start- 
ing either from a fermionico or bosonicEil field theory, 
fictitious Cherns-Simons flux tubes are attached to the 
particles. Comparing these theories, based on singular 
gauge transformations, to our decomposition of the vec- 
tor potential into type-I and type-II components, it is 
obvious that the former theories only take into account 
the singular type-II piece. As a consequence, these mean- 
field theories are not able alone to account for the radial 
component of the Laughlin prefactor. The full Laughlin 
prefactor only emerges when fluctuations of the fictitious 
gauge field beyond the mean field level are taken into ac- 
count. In contrast, our construction correctly accounts 
for fluctuations and produces the full Laughlin wave func- 
tion. In this sense it goes beyond the mean field theorw- 

Another recent theory due to Rajaraman and SondhiEa 
rather artificially attaches complex flux tubes to the elec- 
trons, using a complex Chern-Simons-type vector poten- 
tial The imaginary term in the vector 
potential allows them to obtain the complete a t the 
mean field level. The Aharonov-Bohm phase picked up 
by an electron moving in the presence of such a complex 
vector potential is Qab = —i{^/4>o) J ( a cs —izx a cs )dr. 
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In our construction a moving particle observes real type- 
I and type-II Chern-Simons vector potentials, A" = 
— A l cs — & cs . The type-I component gives the modifica- 
tion = (2w/(po) J a. cs dr while the type-II component 
gives 3?<2 = (2tt/4>o) /(z x a cs )dr. Since !RQ + i^Q = 
—QaBi the two approaches give the same mathemati- 
cal wave function. However the physics is different: our 
approach does not attach bare flux tubes to the parti- 
cles. The full wavefunction is constructed by assigning 
a different role to the two components of the vector po- 
tential. We use the Chern-Simons vector potential a cs 
only as an auxiliary tool in deriving new solutions of the 
Schrodinger equation. The term a cs cancels out from the 
Hamiltonian, in which the total vector potential is A a , 
due only to the uniform applied field. On the other hand, 
the complex vector potential used in Ref. ^3| is an unre- 
alistic characteristic which remains in the Hamiltonian. 



C. Coulomb Tracking 

In previous subsections we analyzed problems with a 
uniform applied magnetic flux density B a . Now we con- 
sider an additional one-body, non-uniform flux density 
.Bi(r), giving the total flux density 



B(r)=B a -B 1 (r). 



(3.22) 



Choosing the associated vector potential in the Coulomb 
gauge, V- A = 0, allows us to interpret B(r) as generated 
by a type-I vector potential. In general, non-uniform per- 
turbations of this kind lift the degeneracy of the LLL, and 
the one-electron eigenfunctions ca nno t be found analyti- 
cally. However, according to Eq. (3.2), if B(r) is tracked 
by a scalar potential Vt(r), i.e. 



HbBi(t) = Vi(t), 



(3.23) 



the degeneracy of the LLL is preserved. Since any 
nonzero spatial average components of the flux density 
and scalar potential only have a trivial level shifting ef- 
fect on the LL structure, we assume that the mean values 
of B\ and V\ are zero in the following treatment. Al- 
though the LLL energy does not change, the basis func- 
tions spanning this level pick up an extra prefactor 



(3.24) 



where SJQi(r) is the solution of the real part of Eq. (3.1) 
with A 1 the solution of V x A 1 = — -Bi(r). 

This result for the single electron states was pointed 
out in Ref. [l7| although all explicit analysis was based 
on a Pauli spinor Hamiltonian with no scalar potential. 
However, we proceed further and suppose that the per- 
turbing magnetic field is due to a two-body, type-I vector 
potential in a system of N electrons. Our goal is to con- 
struct a Hamiltonian, which involves realistic two-body 
scalar potentials, such as the screened Coulomb interac- 
tion, 



,-\ri-Tj\/i 



< 9c < 1, (3.25) 



whose ground state can be constructed exactly. The in- 
teraction is characterized by an exponential screening pa- 
rameter £ and a linear strength parameter g c , with nomi- 
nal value g c = 1. According to the tracking equation Eq. 
( |3.2[ ) this requires the presence of a screened magnetic 
monopolc-like field 



5me 2 e -|r t -r,|/£ 



(3.26) 



with strength g m — g c . With this proviso the degen- 
eracy of the LLL remains intact and the ground state 
wavefunctions have the form 



tf'(r) = e- KQl *(r), 



(3.27) 



wherejftQi({r}) is again the solution of the real part of 
Eq. (O) with V, x A 1 = -B x (vi), and * is the wave 
function in t he ab sence of V\. In the present case, with 
B\ given by (3.26), one obtains 



Ei 



-In 



(3.28) 



where Ei(r) is the exponential integral function Ei(r) = 
— S°°r dtexp(— t)/t, and a = 2m*e 2 g c /h 2 . In the strict 



(3.29) 



Coulomb limit £ — > oo, Eq. (|3.2q) reduces to 



r,: - r, 



j>i 



Assuming tacitly that the ground state structure does 
not change drastically in the presence of a finite electron- 
electron interaction, i.e., the system remains in the 
Laughlin liquid phase and does not cross over into a 
Wigner solid J2J we pick up the Laughlin wavefunction 
from the set of degenerate wavefunctions of the LLL, and 
set ^ = ^ l in Eq. (3.27). In the pure Coulomb case, us- 
ing Eq. ( 3.29| ), this leads to the Coulomb tracking wave 
function 

N N I 

* c =nn - ex p «ee i r » - r i\ 

i j>i \ i j>i 

-jC 2 En 2 )- ( 3 - 3 °) 



The magnetic flux density in Eq. (3.26) may appear 
rather artificial, but it can be viewed as a smeared out 
analog of the standard Chern-Simons flux tube. In 
comparison with the Laughlin wavefunction the extra 
Jastrow factor exp(— 3fQi) pushes the particles further 
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apart, decreasing even more the probability of two par- 
ticles approaching each other, and thus decreasing the 
Coulomb energy. However, this factor does not modify 
the effective filling factor: since appears in the ex- 
ponent, it does not change the polynomial prefactor. In 
a realistic situation with a Coulomb interaction between 
the particles, \&c may be a usefu l var iational ground state 
wavefunction. If one uses Eq. ( 3.2§| ), the corresponding 
screened Coulomb tracking function depends on two pa- 
rameters £ and g c , which can be adjusted variationally, 
independently of their nominal values in Eq. ( 3.25| ). 




FIG. 1. For the model parameter space shown, g c and 
g m are the strengths of the Coulomb and tracking magnetic 
field interactions, respectively. Exact tracking solutions in 
the space of two-body couplings exist along the tracking line, 
g m = g c - Laughlin's wavefunction ^ l is an exact solution at 
gm — g c = 0. The real physical problem is characterized by 
g c = 1, and g m = 0. Any point on the tracking line could be 
used as a variational trial state. 

Figure |l| illustrates schematically the situation in the 
space of model parameters: Exact ground state wave- 
functions can be constructed along the straight line de- 
noting models with equal strengths of the Coulomb inter- 
action and its tracking two-body magnetic field. Laugh- 
lin's wavefunction ^j, refers to the special point g m — 
g c = 0, while the full Coulomb tracking wavefunction 
refers to the point g m = g c = 1. A priori, *S?l 
and are equally good (or bad) representations of 
the wave function for the the real physical case g c = 1, 
g m = 0. Improvement on the FQHE ground state could 
be achieved by considering tracking solutions where the 
common value of the coupling constants g m = g c = g is 
used as a variational parameter, which refers to moving 
along the tracking line of the figure. 

The prefactor of ^c, including the real factor 



exp(— 3?Qi), is not completely anti-analytic; thus it is 
not entirely in the lowest Landau level. Nevertheless, it 
has presumably an extensive overlap with the LLL due 
to the dominating anti-analytic polynomial factor. In a 
situation where the Coulomb energy is comparable to the 
magnetic energy, Landau level mixing is not negligible, 
and $c could be tested directly as a variational ground 
state. Note that the form of t he Ja stro w cor rection fac- 
tor exp(— SRQi) taken from Eq. (3.28) or (3^22) is different 
from the form used recently by Price et alolt would be 
interesting to calculate numerically the variational upper 
bound on the ground state energy set by the tracking 
solution and compare it directly to that of Ref. |24|. 



IV. CONCLUSIONS 

In summary, we analyzed a system of N electrons in- 
teracting with externally applied non-uniform scalar and 
vector potentials, and with each other through two-body 
potentials. Our approach was based on transforming 
the Schrodinger equation into a nonlinear Riccati equa- 
tion and then linearizing the latter with an appropriate 
ansatz. We showed that when the magnetic flux density 
tracks the spatial dependence of a scalar potential, ex- 
act analytic solutions of the Schrodinger equation can be 
obtained in any dimension. These "tracking" solutions, 
which form a Jastrow product, are characterized by a unit 
vector n, and a decomposition of the vector potential. 
In the 2D examples analyzed here, the tracking solutions 
constitute the degenerate ground states of the system. In 
particular, the tracking solution corresponding to the 2D 
electron gas in a homogeneous transverse magnetic field 
was found to be the Laughlin wavefunction whenever the 
vector potential decomposition is based on the Chern- 
Simons vector potential a cs . Our construction, however, 
does not attach singular flux-tubes to the particles, since 
a cs does not appear in the Hamiltonian. 

Adding a realistic repulsive two-body potential, such 
as the screened Coulomb potential, a tracking solution 
exists only if there is a two-body magnetic field present. 
Nevertheless, the resulting tracking wave function may 
be useful in a variational calculation to provide further 
insight into the nature of the fractional quantum Hall 
ground state in the presence of strong electron-electron 
interactions and Landau level mixing. 

In three dimensions there are many possible tracking 
solutions not considered here. Furthermore, the track- 
ing constructions that linearize the Riccati eigenequation 
also linearize the time-dependent Riccati equation, cre- 
ating another domain of exactly solvable problems. We 
leave these interesting problems to future study. 
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